Quantum backaction of optical observations on Bose— Einstein condensates 
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Impressive pictures of moving Bose-Einstein condensates have been taken using phase-contrast 
imaging [M. R. Andrews et ai, Science 273, 84 (1996)]. We calculate the quantum backaction of this 
measurement technique, assuming the absence of residual absorption. We find that the condensate 
gets gradually depleted at a universal rate that is proportional to the light intensity and to the 
inverse cube of the optical wave length. The fewer atoms are condensed the higher is the required 
intensity to see a picture, and, consequently, the higher is the induced backaction. To describe the 
quantum physics of phase-contrast imaging we put forward a new approach to quantum-optical 
propagation. We develop an effective field theory of paraxial optics in a fully quantized atomic 
medium. 
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I. INTRODUCTION 

Bose-Einstein condensates are macroscopic quantum 
states of many atoms with nearly identical wave func- 
tions. In addition to their prospects as laser-like sources 
of matter, the condensates will offer a new fascinating 
testing ground for the fundamentals of quantum me- 
chanics, because of their macroscopic yet quantum na- 
ture. In all experiments performed so far since the pio- 
neering breakthrough |J Bose condensates are observed 
optically. One particularly intriguing method is phase- 
contrast imaging (111. This technique allows to take 
impressive in situ pictures of moving droplets of Bose- 
condensed matter, similar to phase-contrast images of 
delicate living cells 0]. Figure 1 shows a simplified 
scheme of the method. Non-resonant laser light illu- 
minates the sample, travels through, and attains phase 
shifts that are proportional to the density of the conden- 
sate. The remaining part of the apparatus serves to mea- 
sure the acquired phase gradient using the non-scattered 
part of the incident light as a reference. Since the light 
is non-resonant, the condensate is hardly disturbed and, 
indeed, many snapshots or even entire video sequences of 
an individual sample have been taken without significant 
effect g. 

One might suspect g], however, that the observation 
will nevertheless cause a quantum backaction on the con- 
densate. Imagine that phase-contrast imaging were an 
ideal von-Neumann measurement of atomic positions. 
After each measurement the atoms were frozen in po- 
sition eigenstates and their future motion were drasti- 
cally altered. Why is the observed backaction much less 
dramatic? How does phase-contrast imaging affect the 
condensate? The authors of the method Q state that "al- 
though dispersive scattering does not heat up the cloud 
and destroy the condensate, it will change the phase as 
a result of frequency shifts by the ac Stark effect." This 



would still allow "a nonperturbative measurement of the 
number of condensed atoms, which is the variable com- 
plementary to the phase (so-called quantum nondemoli- 
tion measurement)." 

In this paper we find that two processes contribute 
to the backaction. One is the phase diffusion mentioned 
above ||] and the second is a depletion of the condensate. 
The authors of Ref. 0] have already seen indications of 
a gradual depletion. However, they have attributed this 
to a residual light absorption. We find it likely that their 
qualitative experimental findings do in reality indicate 
the effect of the quantum backaction. 

The physical reason for the quantum backaction is the 
local interaction of the meter (the incident light) with 
the object (the condensate). Therefore, to understand 
the quantum nature of phase-contrast imaging, we must 
study the quantum propagation of light in Bose-Einstein 
condensates. To a remarkable degree of accuracy || Bose 
condensates act simply as dielectric media on the inci- 
dent light. In the limit of a large detuning from optical 
resonances of the atoms we can neglect the imaginary 
part of the refraction index, i.e. we consider the con- 
densate a lossless and dispersionless dielectric medium. 
Quantum-optical propagation in bulk, lossless media has 
been studied already in detail q|. A bulk medium, how- 
ever, is hardly affected by light traveling through, and 
the quantum backaction was rightfully ignored so far || . 
In contrast to bulk matter, a delicate quantum gas may 
feel the disturbances caused by incident light. To un- 
derstand actio and reactio of light in quantum gases we 
essay a quantum theory of paraxial light propagation in 
fully quantized, bosonic media. We start from the safe 
ground of canonical electromagnetism in dielectric me- 
dia. Then we develop an approximate theory that allows 
an analytic calculation of the backaction. Our theory is 
inspired by the traditional paraxial approximation |7J] to 
study classical light propagation. We quantize this model 
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and conclude Section II with a discussion of the parax- 
ial field commutator, a critical quantity that enters the 
backaction rates. Section II thus develops a general, phe- 
nomenological theory of non-resonant light propagation 
in matter waves. 

Section III is more specific. Here we derive the master 
equation that describes the backaction of the propagating 
light on the condensate. The light plays a double role: 
It acts both as a meter and as a reservoir. We apply the 
Born-Markov approximation of reservoir theory 8 1 to de- 
rive a general master equation that describes the evolu- 
tion of all (condensed and non-condensed) atoms. Then 
we restrict our attention to the condensate and analyze 
the two dissipative processes that turn out to appear due 
to the quantum backaction: phase diffusion and atomic 
depletion. We conclude Section III estimating the order 
of magnitude of the two effects. We find that the deple- 
tion far outweighs the phase diffusion. Furthermore, our 
estimation seems to indicate that the depletion is experi- 
mentally significant. Condensates with few atoms are es- 
pecially vulnerable to optical observations, in agreement 
with experimental evidence H. 



II. LIGHT AND MATTER WAVES 

A. Canonical theory 

Let us start from the canonical theory of electromag- 
netism in a dielectric medium made of a single matter 
wave ip. The total Lagrangian density of light and mat- 
ter is (in SI units) 
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The first term describes a free electromagnetic field char- 
acterized by the potentials U and A that constitute the 
field strengths 



E 
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-VU, B = VxA. 
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The second term plays a double role. The term describes 
both the effect of the medium \ip\ 2 on the light E 2 and 
the backaction of the light on the medium. The third 
term 
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Lagrangian density j£f describes the physics of a single 
matter wave that interacts non-resonantly with a classi- 
cal electromagnetic field. 

Traditional quantum optics operates mostly with 
Hamiltonians. To find the Hamiltonian density Ji? of our 
model we follow the canonical procedure. We calculate 
the functional derivatives 
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and get 
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DA + — (ip*ip — tp*tp) — J£? = + (5) 



with the atomic Hamiltonian density 
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and the electromagnetic part 
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Observing the Maxwell equation VD = we ignore 
DVC/ = V(DCI) in J#f and arrive at the electromag- 
netic Hamiltonian density 
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This expression is rather unpleasant, since the interaction 
between light and matter appears in the denominator 
(1 + XolV'l 2 )- Of course, we could represent (1 + Xo|V'| 2 )~ 1 
as the geometric series Y^=a(~ Xo|V'| 2 )' y i provided that 
the series converges. The lowest-order theory 
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has a familiar appearance in quantum optics, because 
the Hamiltonian density J^r contains the biquadratic in- 
teraction term D 2 |i/>| 2 . However, if we were to use the 
lowest-order J^f R in Section III to calculate the quantum 
backaction in phase-contrast imaging, we were to face an 
infinite result that clearly contradicts experimental evi- 
dence Can we formulate light propagation in matter 
waves in a different way? 



characterizes the center-of-mass motion of the matter 
wave in the external potential V. One can show that the 
Euler-Lagrange equations generated by minimizing the 
action J _S? d 4 x are indeed Maxwell's equations of light 
in a dielectric medium with susceptibility eoXo|V-'| 2 - Ad- 
ditionally, we obtain the Schrodinger equation of a mat- 
ter wave moving in both the external potential V and in 
the optical potential — eoXoE 2 /2. This proves that the 



B. Optical Schrodinger equation 

Our canonical field theory of light and matter contains 
Maxwell's equations in the form of the Lagrangian _§f 
or, equivalently, in the Hamiltonian Jif. In classical op- 
tics, however, one hardly uses Maxwell's equation to de- 
scribe paraxial light propagation, but applies a charac- 
teristic approximation which is an optical analogue of the 
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Schrodinger equation. Let us sketch a brief derivation of 
the optical Schrodinger equation Q). We start from the 
wave equation 



V^E- 
VE = 
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(a familiar consequence of Maxwell's equations). We as- 
sume the medium Xo IV' 1 2 to be weak and to be gradually 
varying in space. (A mirror, for example, is excluded 
from this model.) We neglect the polarization mixing 
V(VE) and consider paraxial propagation in the z direc- 
tion with carrier frequency ujq. In the following we will 
often refer to the wave number k — loq/c and to the 
wave length A = 27r/fco- We use the ansatz 



E^ -1- ) = £ exp(ikoz — iujot) 
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for the positive frequency component E^ + ) of the electric 
field, and differentiate 
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neglecting the slow variation of the envelope £. Further- 
more, we regard 
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with being the transversal Laplacian d 2 /dx 2 + 

d 2 /dy 2 . Finally, we utilize that (d/dz+d/dct) exp(ikoz — 
iujot) vanishes and arrive at the optical Schrodinger equa- 
tion 
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The complex function tp represents the optical field in 
units of the vacuum noise (Kwo/2eo) 1 / 2 . One can easily 
verify that the Lagrangian density 
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generates the optical Schrodinger equation (|15|) as an 
Euler-Lagrange equation for tp. This proves that ££f 
is a valid Lagrangian to describe paraxial optical propa- 
gation. Lagrangians are defined up to a positive pref- 
actor (because we are only interested in their mini- 
mum). We have chosen the prefactor of ££f such that 
the coupling term ^oXolVfM 2 / 2 = e Xol^| 2 |^| 2 ~ 
eoXo(^ + ' + E^') 2 \tp\ 2 /2 is identical to the interaction 
term eoXoE 2 ^ 2 /2 in the canonical Lagrangian density 
(|l|) . This is necessary, because the term plays the double 
role of describing actio and reactio of light and matter. 

In order to find the Hamiltonian density we follow the 
royal road of canonical field theory. We calculate the 
functional derivatives 
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and obtain 



^fp = — (tp*tp — tp*tp + ip*ip — ip*t/>^ — -Sfp 

= Jzm + 'ftp , 

ihc (dip* f dip\ he 



^0 I |2| / |2 

-— xom m 



(19) 



The last term, —huioXo\f\ 2 \' l { } \ 2 /2, has the familiar ap- 
pearance of a quantum-optical light -matter interaction. 
The other terms of M'f describe the free paraxial propa- 
gation of the optical field. 



More details about the optical Schrodinger equation (and 
beyond) are elaborated in the comprehensive paper M by 
Marte and Stenholm. 



C. Canonical theory of paraxial propagation 

Let us describe the paraxial propagation of light in a 
canonical field theory. We consider a fixed polarization, 
i.e. a scalar electrical field strength E^ + \ and express 
£(+) as 



D. Quantum statistics 

So far we have formulated the canonical theory of a 
single matter wave that interacts non-resonantly with 
an optical field. Let us turn to the description of a Bose 
gas of many atoms interacting with photons. The many- 
body problem is conveniently formulated in the language 
of Second Quantization where operators ip and tp describe 
the matter and light field, respectively. The atoms are 
bosons 
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(20) 
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To find the Hamiltonian H of our model we integrate the 
Hamiltonian density Jfp over three-dimensional space 
and apply partial integration. We obtain 



H = H R 
Hr = Tic J <p 
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The Hamiltonian Hr describes free paraxial light propa- 
gation, Hi accounts for the interaction of light and mat- 
ter and Hm describes the quantum gas of the atoms. 
We include atomic collisions in our model by adding an 
atom-atom interaction term to the potential V, and ob- 
tain for example [Bl 



Hm = 



V> f V 2 V> + Vtp^ + - tp^ip 2 

2m 2 
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Note that our results turn out to be independent on the 
explicit form of Hm, as long as Hm is able to maintain 
a Bose-Einstein condensate. 

An important issue in our theory is the optical field 
commutator 



C(xi -x 2 ) = [^(x 1 ),^ t (x 2 )] 
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If C(x) is a three-dimensional delta function we obtain 
the optical Schrodinger equation ( |l5| ) for ip from the 
Heisenberg equation 
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Is the commutator a delta function? The field operator ip 
is one polarization component of the photon absorption 
operator V(x, t) of Mandel and Wolf jl0|]. The commu- 
tator C(x) is then given in terms of the density <i/i(k) of 
the employed optical modes as Q 



C(x) 
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The crucial point is that the mode density d/i(k) is re- 
stricted, since we consider only those modes that obey 
paraxial propagation. Therefore, we cannot regard C(x) 
as a three-dimensional delta function, in general. How- 
ever, the commutator behaves like 5^ 3 ' (xi — x 2 ) on parax- 
ial test functions that are supported on the spectrum of 
C (that share the same region in k space). To see this we 
introduce the Fourier transformations C(k) and /(k) of 
the commutator and a test function /(x), and calculate 



/c (x , 
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This is the defining property of a three-dimensional delta 
function. In particular, the field operator (p itself is a 
paraxial test function, and thus we can safely regard C(x) 
as (5' 3 )(xi — x 2 ) in integrals containing (p. In this way 
we obtain indeed from the Heisenberg equation (^) the 
optical Schrodinger equation for ip, as is easily verified. 
This shows that our quantum theory of paraxial light 
propagation is consistent with the classical theory. 

On the other hand, for test functions with a broad 
spatial spectrum we put 



dn(k) = 6(oj(k)/uJo - l)d 3 k 

= 5{k z /k Q + \{kl + k 2 y )/kl-l) d 3 k. 

We obtain, for example, the commutator per se 
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see Appendix A. The commutator resembles the propa- 
gator of a free-particle wave, because C satisfies 
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Since commutators are Green's functions, we would ex- 
pect this property from a paraxial quantum optics. 

We have thus seen that the commutator acts differently 
on paraxial and on broad-spectrum test functions. The 
two faces of the commutator will play two different roles 
in the quantum backaction of phase-contrast imaging. 



III. QUANTUM BACKACTION 

A. Phase— contrast imaging 

Phase-contrast imaging is a method to see transpar- 
ent objects (living cells, for example) that are otherwise 
invisible. A transparent body does not absorb incident 
light but imprints a phase shift that is proportional to 
the local density of the object. Manipulations in the focal 
plane of the observing lens system transform the image 
of a phase object into an equivalent image of an ampli- 
tude object Q, see Fig. 1. In the absence of absorption 
the light will not significantly affect a macroscopic body, 
but it might affect a delicate quantum gas. 
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Formulated in the language of quantum measurement 
theory the light is a meter of the atomic object: The 
incident light interacts with the atoms, becomes entan- 
gled with the atomic sample, and is finally measured in 
a classical apparatus (lens system and detectors). The 
light field is also a reservoir: The incident beam, say a 
plane wave, is scattered on the object into many emerg- 
ing modes that form an image. The degree of collecting 
the scattered modes and extracting the imprinted phase 
information determines the performance of the measure- 
ment. The overall performance is ultimately limited by 
the quantum-optical phase fluctuations of the incident 
light and by the performance of the apparatus. 

On the other hand, the average quantum backaction 
(averaged over all experimental runs) is solely determined 
by the entanglement between object and meter. This 
means that the average backaction is entirely indepen- 
dent on the performance of the measurement. We can 
average over the meter to calculate the backaction on the 
object. This averaging leads to a dissipative dynamics 
of the object, formulated as a quantum master equation 



B. Master equation 

We describe the state of all atoms by the density ma- 
trix pi in the interaction picture 



Pi{t) 



jH M t/h 



p(t) e 



(32) 



and employ standard reservoir theory JsJ to derive an 
equation of motion for pi. The interaction between light 
and matter is weak and, consequently, we can apply 
second-order perturbation theory (the Born approxima- 
tion Q). The incident light travels through the object in 
a time t c that is simply given by the extension of the con- 
densate divided by the speed of light. After r c we do not 
expect the light to interact further with the atomic sam- 
ple. Therefore, the effective memory of the light-matter 
interaction is roughly given by r c and is very short. So, 
in addition to the Born approximation, we can apply the 
short-memory (Markov) approximation using an integra- 
tion time To that is much larger than the time of flight 
t c but is still smaller than a characteristic time of the 
atomic sample. We will see later in Sections C and D 
that our results do not depend on the precise value of 
the memory time To, provided that To exceeds r c . This 
justifies a posteriori our Markovian theory. We obtain in 
the Born-Markov approximation fefl 
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~dt~ ~ h 
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with 
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tT R [H n (t), [Hii{t'),pi{t)p m (t)}} dt' (33) 
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where ip and <p are Schrodinger-picture operators, and 
Pir denotes the density matrix of the light field in the 
interaction picture 



Pin{t) = e 



,-iH R t/h 



(35) 
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The incident light is a plane wave 

tr R {(^ t (xi)^(x 2 )/5 i? } = — — l. 

tx R {^> t (xi)^ t (x2)^(xi)^(x 2 )pfl} = const, (36) 

that travels in the z direction and carries an intensity 
(energy flux) of /. Apart from the value of the intensity 
we do not need to specify the quantum state of the illu- 
minating plane wave. The atomic sample is a thin phase 
object, thin enough to neglect the optical diffraction in- 
side (the V^_/2fco term in Eq. (22) of the Hamiltonian 
H R ). This means that (pi obeys the equation 



d 1 d 
dz + cdl) Vl 







with the obvious solution 

ipi(x,t + to) = ipi(x- cte z ,t ) . 
The total number of atoms is fixed, i.e. 

[p, { tP^d 3 x] = 0. 



(37) 



(38) 



(39) 



We assume the existence of a Bose-Einstein conden- 
sate with wave function tpo(x.,t). To describe this phe- 
nomenon, we expand the atomic annihilation operator 
ipi{yL,t) in the interaction picture into a complete, or- 
thonormal set of atomic modes 



i>i{^t) = y^„(x, t)a v 



■0*(x,t)'(/; ly (x,t)d 3 a; = 5^ 



5>i(*i, t)^(x 2 ,<) = ^( Xl -x 2 ) 



(40) 



We assume that only the ipo mode (the condensate) is 
significantly populated, d v pi = for v ^ 0. This means 
that we consider the case of zero temperature and that 
we neglect the fluctuations Q of the above-condensate 
part, for simplicity. For repulsive atom-atom interac- 
tions the spatial wave function ipo of the condensate is 
dominated by the balance between the interatomic re- 
pulsion and the expernal potential H. The state of the 
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condensate is described by the reduced density matrix 
Po = t^ACPl- 

For finding the master equation of the condensate state 
Po, we average the master equation (|33| ) with respect 
to the above-condensate part. This procedure requires 
some lengthy yet straightforward calculations that we 
prefer to present in Appendix B. We obtain the result 



—7- = -teAC^Pi = -Apo - C2P0 , 
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(41) 
(42) 
(43) 
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V 

T L = J G(cre z ) dr . 
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(44) 
(45) 

(46) 



We see that no macroscopic light force acts on the con- 
densate. This is easy to understand, because light forces 
are caused by intensity gradients, and in our case the 
incident light is a plane wave with uniform intensity /. 

We notice that two dissipative processes, Cipo an d 
C-2,po-, contribute to the backaction. The first one is a 
phase diffusion (see Section C) and the second one is a 
bosonic depletion (see Section D). What is the physical 
origin of the two processes? 

For forming an image of the condensate, the incident 
light is scattered on the atoms. Therefore the light trans- 
fers momentum (and energy) to the atoms. After the 
scattering a condensate atom may find itself in the con- 
densate again or it may be expelled, with the (complex) 
rates Tp and Tl — Tp, respectively. The atom that re- 
joined the condensate will suffer from a loss in phase co- 
herence, causing phase diffusion, and the expelled atoms 
will gradually deplete the condensate. Let us look into 
the two processes in more detail. 



C. Phase diffusion 

Let us study the phase-diffusion Cipo in the atomic 
number-state basis of the condensate 



(m\Cip \n) = [ilmTp (m 2 - n 2 ) + Rer P (m - n) 2 ] 



x (m\p \n) 



(47) 



We see that the imaginary part of Tp governs a Hamil- 
tonian process, i.e. a commutator of po with a Hermitian 
effective Hamiltonian that is proportional to the atom 
number. We may call it a Kerr-type process. 



Let us turn to the complementary picture where we 
describe C\po in the basis of phase states |l5| \4>) = 
(2-7r) _1 Y^=o ex P(* n 0) \ n )- We obtain 
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The first term is clearly a diffusion process of the atomic 
phase (and hence justifies our terminology of calling C\po 
a phase diffusion). The diffusion rate 7^ is the real part 

of r P . 

To calculate jp we utilize that the characteristic time 
of the condensate is much larger than the time of flight 
r c of the incident light and the integration time To- This 
means that we can regard \ipo(x, t')\ 2 as being equal to 
|V>o( x ) t)\ 2 m the expression (E3) of tp, i.e 



|^o(x,i')| 2 = l^o(x,i)| 2 =Po(x,i). 
We apply the Fourier transformations 
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p (k)^ J Po ( X )e-^d 3 x : 



G(k) = / G(x)e- Jkx d 3 x , 
and use the convolution theorem to obtain 
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(50) 



1p = 



2(2tt)= 



|p (k,t)rG(k)e lcr ^^fcdr. (51) 



Because cto is large compared with the extension of the 
condensate we can approximate 
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|po(k,i)| 2 G(k),5(fc z /fc )d 3 fc. (52) 



We utilize the definition (|46|) of the G function and find 
that 



2 

G= ^^/G(k + fc e 2 
4 nc 



(53) 



with G(k) being the Fourier-transformed commutator 
(^H) expressed in terms ( p7|) of the mode density d/i(k). 
We obtain 



IP = 



(54) 



Since a typical Bose-Einstein condensate extends over 
many optical wave lengths, the Fourier-transformed spa- 
tial distribution po(k) is narrow compared with the op- 
tical mode density. Furthermore, in Eq. (|54|), j?o(k) is 
moved by koe z to an area in k space where it overlaps 
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with the optical field. The spatial distribution of the con- 
densate thus probes the paraxial optical mode density. 
In this case, see Section IID, the optical-field commuta- 
tor acts as a three-dimensional delta function. We set 
d/i(k) — d 3 k and obtain, finally, 
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\po(k x ,k y ,0,t)\ 2 dk x dk y , (55) 
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7p = — —1a 




Po(x-,t)dz dxdy. (56) 



In this way we have expressed the phase-diffusion rate "fp 
in terms of the spatial probability distribution po (x, t) for 
a single condensed atom. 



D. Depletion 

The second backaction process of phase-contrast imag- 
ing is a familiar bosonic depiction. The Liouvillian C2P0 
contains a Hamiltonian phase shifting that is propor- 
tional to ImTx — ImlV and a depletion that occurs at 
the rate jl — 7p with jl = Rel^. We have calculated 
7p in the previous section, let us turn to jl here. First, 
we proceed along similar lines as in the calculation of "fp 
and get 



1 



„ = Z[X0r 

1L A he (2tt) 3 



6(k z /k - 1) d/i(k) . 



(57) 



Because 



|Po(k) 



p (x)e 



lkx d 3 a 



we see from Eqs. (p3) and (p7]) that 



Jl > 7p ■ 



< / p (x)d 3 x = 1 (58) 



(59) 



The overall depletion rate jl — 7p is non-negative. This 
satisfactory result illustrates the consistency of our the- 
ory. 

Note, that we cannot replace the mode density dfi(k) 
in Eq. ([[t]) by d 3 k, as we have done in the case of phase 
diffusion, because in Eq. ( |57| ) we integrate d/x(k) over a 
broad spectral range. Therefore, we must use the broad- 
spectrum behavior (^o|) of the commutator (^7|) to cal- 
culate jl- We insert the definition ( ff6| ) of G into the 
expression ( filf ) of Tl and obtain for 7^ = Rel^ the in- 
tegral 



1L 



ko Xoj 
8 he 



+ CT 



C(ze z )e 



—ik,QZ 



dz . 



(60) 



Then we apply the explicit formula (30) of the commu- 
tator and get 



1L 



+ CT 

— X0j\-3 f dZ 

Aihc J z — iO 

-CT 

— ^I\- 3 — ( I — 



2 he 



2ni 



z-iO 




(61) 



with the contours indicated in Fig. 2. We use the residue 
theorem for the first integral and calculate the second 
explicitly. We obtain, finally, 



2 2 



(62) 



Our result is independent on the value of the integra- 
tion time ro and is therefore consistent with the assumed 
Markovian behavior of the quantum backaction. Further- 
more, Eq. ( |62| ) does not depend on the spatial shape of 
the condensate. The loss rate is universal. 

Equation (|6(]) shows that the depletion rate 7^ de- 
pends on the optical-field commutator C(x,y,z) at 
x, y = integrated over z. So, if the commutator were 
a three- dimensional delta function, as in usual quantum 
optics, our calculations would give an infinite backaction, 
a result in clear contradiction to experimental evidence 
0. At this point the reader will appreciate the reason 
why we have developed a quantum theory of paraxial 
propagation in Section II. This theory employs both a 
biquadratic interactio n d23|) of light and matter and a 
modified commutator (|30|)that produces a finite backac- 
tion. 

On the other hand, a first-order perturbation theory 
(^) of canonical electromagnetism (||) is also biquadratic 
in the light-matter interaction. Here the canonically con- 
jugate variables are A and D and hence we should postu- 
late that the commutator of A and D is proportional to 
a transversal delta function. However, this approach pro- 
duces an infinite backaction. Why is perturbation theory 
not appropriate in the canonical Hamiltonian (||)? Re- 
member that here we approximate the light-matter in- 
teraction (l + XolV'l 2 ) -1 by the first two terms of the 
geometric series 2Z^Lo( — XolV^I 2 )^- If "0 is an operator 
tfj, as it should be in the quantum statistics of matter 
waves, the particle nature of matter will make the den- 
sity operator potentially large, regardless how small 
Xo is. (For example, the density of a localized particle is 
a three-dimensional delta function.) Therefore, we must 
not approximate (l + xoi> 'VO -1 by 1 — Xo^^ an dj more- 
over, we must not see (1 +Xoi } ^/O -1 as a geometric series 
at all. 

Let us add a more intuitive picture to our mathe- 
matical arguments. Phase-contrast imaging involves a 
particle-particle scattering of atoms and photons that 
depends on the ability of the incident photons to lo- 
calize. The optical-field commutator describes potential 
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photon fluctuations and therefore photon localization as 
well Jl6|| . Mandel and Wolf [[l6| discuss in great detail 
that photons can manifest themselves as particles only 
within a volume that is bounded by the optical wave 
length. The finite localization volume of photons thus 
leads to a finite quantum backaction. 



E. Estimation of the effect 

How large is the backaction? Let us estimate solely the 
order of magnitude. First, we compare the depletion rate 
7l with the phase-diffusion rate jp. The latter depends 
on the spatial shape of the condensate. To find a rough 
estimate for 7p we assume that the condensate extends 
over a x , a y , and a z in x, y, and z direction, respectively, 
and we model the Fourier-transformed probability dis- 
tribution by the Gaussian 



p (k) = exp {-\a 2 x k 2 x - \a 2 y k 2 y - \a\k 



We obtain from Eq. (B5[) the phase-diffusion rate 



IP 



4 he \{2'Ka x ){2'Kay) 



(64) 



or, expressed in terms of the depletion rate 7l via Eq. 



A 2 



7P = 7L 



(2ira x )(2ira y ) ' 



(65) 



A typical Bose-Einstein condensate extends over many 
optical wave lengths A. Therefore, 



1L > 7p 



(66) 



Depletion far outweighs phase diffusion in the backaction 
of phase-contrast imaging. 

The depletion rate (p2) is proportional to the intensity 
of the illumination. A moderate intensity will thus lead 
to a low backaction. On the other hand, low intensity 
light exhibits large quantum-optical phase fluctuations. 
To give a rough estimation, we use the uncertainty prod- 
uct Wh ' 



)fa K 1 



(67) 



of the phase and photon-number fluctuations. A 
coherent-state illumination has Poissonian photon statis- 
tics with a variance S 2 n that equals the mean n, i.e. 



5 2 <f> « n" 1 . 



(68) 



The mean photon number is the product of the average 
photon density I/(Tlloqc) with the optical mode volume 
during the observation time At. The mode volume is 
roughly given by a cylinder of radius A and length cAt. 
Therefore, we obtain 



/ o . X 2 IAt 
— ir\ cAt = ir— 



(69) 



One can see an image of a phase object only when the 
produced phase shift A(f> exceeds the phase noise 8(f). Let 
us estimate the signal A<f> for a condensate of N atoms 
with single-atom probability distribution pq (x) : 

A<j) = J (k - k )dx 

Vl + XoiVpo(x)-l) dz 
-XoAr; (70) 
where r\ denotes the effective two-dimensional density 
■q = / p (x)dz . (71) 



(72) 



2tt 

T 

7T 

A 



(63) 

For our simple model (|63l) the largest value of ?] is 



77 = 



2T7a x ciy 



Let us compare the signal-to-noise ratio A(j)/5(f> with 
the induced backaction. According to our theory, the 
incident light depletes gradually the condensate during 
the observation time At. The average number of atoms, 
tr{aQaoPo} 5 decays by the factor of exp(— 2^LAt). We 
express our result (|6^) for the backaction rate 7l in terms 
of the signal-to-noise ratio A4>/84> using the estimations 
and find 



(73) 



The more atoms are condensed the lower is the re- 
quired intensity for producing a faithful image and, con- 
sequently, the lower is the induced backaction. On the 
other hand, the more macroscopic the wave function is 
the lower is the probability density 77 and the larger is the 
backaction rate. We also realize from Eq. ([73|) that the 
depletion factor exp(— k) is independent of the suscepti- 
bility xo (and in particular of the detuning from atomic 
resonances), given a fixed signal-to-noise ratio A<f>/5<j>. 
This is easy to understand, because for a l arg e xo the 
signal Acj) is large, but so is the backaction (p2\). 

We obtain for the critical value A<p w 5(f) the deple- 
tion constant k w (NX 2 rj)~ 2 . Let us assume that a x a y 
is roughly 10 4 A 2 . In this case the depletion constant k 
is about A~ 2 10 10 . A small condensate [|| with N « 10 3 
atoms disappears immediately! Therefore, a continuous 
monitoring of the condensate is impossible, in agreement 
with experimental observation j|. Larger condensates 
are more robust. For N sw 10 6 atoms || wc obtain a 
decay factor exp(— n) of about 1 — k w 1 — 1CP 2 . The 
depletion is in the order of one percent which agrees well 
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with the experimental observation [0. However, the au- 
thors of Ref. have interpreted the decay of the con- 
densate as an residual absorption effect. Our calculations 
seem to indicate that the quantum backaction of phase- 
contrast imaging has caused the observed depletion via 
the momentum transfer of the illuminating light quanta. 

Finally, we remark that the use of appropriately 
squeezed fl§| illumination would enhance the signal-to- 
noise ratio at a given, moderate quantum backaction. 



IV. SUMMARY 

Phase-contrast imaging induces a quantum backac- 
tion on Bose-Einstein condensates. The backaction con- 
sists of two dissipative processes: a phase diffusion and 
a gradual depletion of the condensate. If phase-contrast 
imaging were a quantum nondemolition measurement of 
the condensate density, solely phase diffusion would oc- 
cur. We have shown, however, that the depletion is the 
dominant process. Our estimations indicate that the 
quantum backaction is indeed an observable effect that 
limits the otherwise destructionless character of phase- 
contrast imaging Condensates with a few number of 
atoms are especially vulnerable to optical observations. 

In our theoretical study we ignored entirely the absorp- 
tion of light (the imaginary part of the refractive index) 
and we focused our attention on the quantum backac- 
tion. This is justified in the limit of strong detuning, be- 
cause here the absorption virtually vanishes. Note that 
our theory is only valid as long as the condensate is still 
significantly populated. We have assumed a priori the 
existence and the dominance of the condensate. 

We developed a quantum theory of paraxial light prop- 
agation in many-atom samples that are off-resonant with 
respect to the light. This theory goes beyond the immedi- 
ate purpose of the present paper and may find wider ap- 
plication in the fascinating subject of combining concepts 
of quantum measurement theory with quantum-optical 
propagation and the statistical theory of quantum gases. 
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APPENDIX A 

In this appendix we calculate explicitly the optical field 
commutator C(x) and we discuss a few general properties 
of C(x). We start from the expression ( p7j) of C(x) in 
terms of the mode density dfj,(k), 



C(x) 



1 



ikx 



dfj,(k) , 



(27T) 3 

dn(k) = 6(k z /k + l(k 2 x + kl)/kl 



1) d 6 k . 



(74) 



To calculate the Fourier integral we introduce dimension- 
less variables, k x = kg (^ x , k y = k C, y with 



k z =k -±k {g + Q, 



(75) 



and obtain 



C(x) = 



I-'-'' 
ft o 

(2tt) 3 



-iQk (z-i0)/2+ik o xC. 



x d( x 



e -iC 2 y k {z-iO)/2+ik oy ( s d Q y ( 76 ) 



In Eq. (JT^) we have given z an infinitesimally small yet 
negative imaginary part —iO to ensure that the integrals 
converge. The calculation of Eq. (76) is straightforward, 
and we obtain the result 



C(x) 



h 2 



i{2n) 2 (z -iO) 



exp 



ik 2 



y 



ikoz 



(77) 



Let us discuss a few elementary properties of the paraxial 
commutator. First, we see from the mode density in Eq. 
([74]) that C(x) obeys the partial differential equation 



V 2 
2fc 



oz 



(78) 



The commutator behaves like a paraxial monochromatic 
wave that travels freely in space without being influenced 
by a medium. Furthermore, when z approaches zero, 
the commutator C'(x,y,z — > 0) vanishes effectively for 
x, y 0, because the exponential function in formula ( |7^ ) 
is highly oscillating in this limit. To see what happens 
at x, y = we integrate C(x, y, z — > 0) with respect to x 
and y. Since 



exp 



2vr 



ik , 2 
2z V 



o 

z 

2m — 

k 



-y 2 ) 

iko . 
ex Pl — r 



dx dy 



we obtain, finally, 

C(x,y,z 



0) 



2n 



S(x)S(y). 



(79) 



(80) 



The commutator C(x) appears as a paraxial Green's 
function. 
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APPENDIX B 

In this appendix we derive the master equation ( pd| ) 
that describes the quantum backaction of phase-contrast 
imaging on Bose-Einstein condensates. We start from 
the general master equation ([53|) and proceed in two 
steps. First, we reformulate the general master equa- 
tion ([33]) utilizing the fact ( pf ) that the incident light 
is a plane wave and using the conservation ( |39|) of the 
total number of atoms. Then we average the resulting 
Liouvillian with respect to the above-condensate part. 
This, finally, gives us the desired master equation for the 
condensate. 

To begin, let us first consider the Hamiltonian term 
tr R [pipi R ,Hn] in Eq. (|33| ) that describes the average 
light force on the atomic sample. We calculate the com- 
mutator of the interaction Hamiltonian and the total den- 
sity matrix in the Schrodinger picture using Eq. (|3^) and 
the conservation of the total number of atoms, Eq. (39), 



tr R [Hi,pp R ] = ^Xo [ / fti>d 3 x,i> = 



(81) 



Consequently, we get in the interaction picture 

ti R [H u ,PiPm]=0 ■ (82) 

The overall light force of an uniform illumination van- 
ishes. 

Let us turn to the Liouvillian part Cpi of the master 
equation (^3|). First we calculate the correlation func- 
tions (|36| ) of the incident light in the interaction picture. 
We neglect optical diffraction inside the atomic sample, 
see Eq. @, and apply Eq. @ to get 

trfij0j(xi,ti)0/(x2,f 2 ) pm(t 2 )\ 

= tr R j<£j(xi - c(ii - t 2 )e z ,t 2 )0i(x 2 ,t 2 ) Pm(t 2 )\ 

= tr^{e^* 2 ^( Xl _ c ( fl - t 2 )e z )<p{ X2 ) PR {t 2 ) e tflRt *} 

= tr R {^(xi - c(ti - t 2 )e z )(p(yL 2 ) PB.(h)} 
I 



exp[ik a (z2 -zi + c(ti - t 2 ))) 



huioc 

and along similar lines 



(83) 



tr fl |(^|(xi, ti)ip\(x.2, < 2 )^/(xi, ti)ipi(x. 2 , t 2 ) /5/_r| 



= const 



(84) 



To evaluate the Liouvillian Cpi in Eq. (|33j) we calcu- 
late the double commutator using the shorthand notation 



2 , ,2 



h ui. 



-Xo / / tv R {$$1 0101 ^2^2 ^ 



2 pIpIR 



~ Wl 0101 PIPIR V#2 0202 

- VW2 0202 PIpIR ipypi 0101 

+ PIPIR i> 2 i>2 0202 ^1^1 0101 1 

x d 3 x\ d 3 x 2 



-1, ,2 



-Xo 



$$1 ^2 pi - ^\i>2 Pi i>\i>\ 



tl'fi {0! 



Ifl ip 2 f 2 PIR 



}d 3 x 



x d 3 x 2 + H.c. 



(85) 



Then we utilize the commutator relation of the radiation 



field fl25| ) and the diffraction-less propagation (38) in the 
interaction zone to write tr R {(p\0i (p\(p 2 Pir} in normal 
order, 



trjjj^Ivi ¥>2<r~2 I'll! 



} 



020102 P/flj 



+ tr fl [ip\02 Atr} C(xi - x 2 - c(<i - t 2 )e z ) . (86) 

The normally-ordered correlation function (Q) is con- 
stant, and so tr R {(p\(p\(p~\(p 2 pi R } does not contribute to 
the double commutator (|85|), when we take into account 
the atom-number conservation ( |39|). To proceed with the 
remaining term, we utilize Eq. (Bq) to express the cor- 
relation function tT R {0\0 2 Pir}, define in Eq. ( |46] ) the 
function G(x), and get 



Cpi = 



Wi ■4>ii>2'>p2Pi ~ WaPzWi 



x G(xi - x 2 - c(i - t')e z ) d 3 xi d 3 x 2 dt' + H.c. . (87) 

In this way we have found the Liouvillian that describes 
the quantum backaction of an uniform illumination on 
the total atomic sample. 

So far we have not made any essential assumption on 
the state of the atoms and in particular we have not ex- 
ploited the fact that the atoms are Bose-condensed. To 
incorporate Bose-Einstein condensation we expand the 
atomic annihilation operator ipi in the interaction pic- 
ture in a complete, orthonormal basis of atomic modes, 
see Eq. ([it)]). We assume that only tpo is significantly oc- 
cupied and average with respect to the above-condensate 
part being in the vacuum state. (We utilize that a v pi 
vanishes for v 7^ 0.) For this we must calculate some 
atomic correlation functions that occur in the Liouvil- 
lian (|87j). We abbreviate tiACPi by po, and find 



tr ac ^l^l^i^ P/| 

= (dt) 2 dl |?/>o(xi,ii)| 2 |-0o(x 2 ,t2)| 2 po 
trAcj^I^iP/j = aJao|V'o(x, t)\ 2 p . 

tl\4c{^2^2 Pi $$1} 



(88) 
(89) 
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= tr AC Upl&o pi al4>%\ V , o(x 2 ,t2)V'o(xi,ii) , (90) 
tr ac |^2«o /O/ aJV'i} 
= y^tr^c |sJ,aopjaJa (tl |^*(x2,t2)V'M( x i'*i) 

= ^ tT AC {al&Q pi ala v \ ^*(x 2 , t 2 )tp u (xl, ii) 

= aJ«o Po 4«o V'o( x 2^2)'0o(xi,ii) 

= a^ao poalaoipQ{x 2 ,t 2 )^o{^i,ti) + a 0/ o aj <5 (3) ( x i ~ x 2) 
-a /5 aj Vo( x 2:*2)V'o(xi,ti) . (91) 

We apply the results (|8^-^l|) to obtain after averaging of 
the Liouvillian ( |S7| ) with respect to the above-condensate 
part the final master equation (pl|). 
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FIGURE CAPTION 



Figure 1: Scheme of phase-contrast imaging |Q. 
To produce an image / of a transparent object the 
non-scattered light is phase-shifted (phase-contrast mi- 
croscopy) or blocked (dark-ground method 0) in the 
focal plane F of the magnifying lens L. 

Figure 2: Integration contours of Eq. ( |6l| ) in the com- 
plex z plane. The closed contour (1) & (2) encircles a 
(z — i0) _1 singularity (dot). 
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